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ABSTRACT 

We demonstrate that the combination of GALEX UV photometry in the FUV (~ 1530 
angstroms) and NUV (~ 2310 angstroms) passbands with optical photometry in the 
standard U, B, V, R, I filters can efficiently break the age-metallicity degeneracy. We 
estimate well-constrained ages, metallicities and their associated errors for 42 GCs in 
M31, and show that the full set of FUV, NUV, U, B, V, R, I photometry produces age 
estimates that are ~ 90 percent more constrained and metallicity estimates that are 
~ 60 percent more constrained than those produced by using optical filters alone. The 
quality of the age constraints is comparable or marginally better than those achieved 
using a large number of spectroscopic indices. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: fundamental parameters 
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Broad-band photometry has traditionally played a key role 
in the study of a variety of astrophysical objects. Often, it 
is of interest to decipher the formation histories of various 
systems, e.g. globular clusters (GCs), by age-dating them. 
The vast majority of stu dies have so fa r been restricted t o 
optical photometry e.g. (lYi et alJliool ; IPeng et al.ll2004h . 
However, the use of optical photometry has been plagued by 
the well-known age-metallictty degeneracy (AMD) whereby 
young, metal-rich stellar populations produce optical colours 
which are indisting uishable from those due to old metal- 
poor populations (I Worthevl Il994) . It is therefore difficult 
to simultaneously pin down the age and metallicity of an 
object using optical colours alone. While combinations of 
metallicity-sen sitive and age -sensitive colours, e.g. U — B 
versus B — V (|Yi et alj 120041 ). can produce rough estimates 
of age and metallicity, the errors on these estimates can be 
extremely large. For example, the age uncertainty cou ld be 
as large as the age of the universe (e.g. lYi et al.1l2004l )! 

Ages can be better studied, and the AMD alleviated 
spectroscopically, by combining a sufficient n umber of age 
and metallicity-sensitive line indices (see e.g. iTrager et al.l 
l200d : IProctor fe Sansomll2002l ; ICaldwell et alj|2003h Forex^ 
ample. iTrager et al.l 1 200ol ) use the Hp, Mg 6 and <Fe> line 
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strengths to estimate well-constrained SSP-equivalent ages 
and metallicities for a sample of local elliptical galaxies. Sim- 
ilarly, IProctor fc Sanso ml (|2002h use a combination of 24 line 
indices to break the AMD and produce SSP-equivalent pa- 
rameters with small uncertainties for the central regions of 
32 systems taken from both the Virgo cluster and the field. 
However, such methods rely on accurate measurements of a 
sufficiently large ensemble of indices for each object of in- 
terest. 

Rest-frame UV flux, shortward of ~ 3000 angstroms, 
provides an attractive alternative. From a theor etical point 
of view, it has already been sug gested (e.g. Iy] |2003|) that 
the UV spectral ranges may provide a better handle on the 
ages of objects and thus reduce the AMD. It is the aim of 
this Letter to quantify the usefulness of the UV in breaking 
the AMD and providing robust age estimation, especially in 
the context of simple populations like GCs. If the addition 
of UV photometry were to provide comparable constraints 
in the age-metallicity parameter space as, for example, large 
numbers of spectral line indices, then the advent of large- 
scale UV photometry from the GALEX mission (Martin et 
al. 2005) could provide an accessible means of age-dating a 
wide variety of nearby objects using photometry alone. For 
example, in its Medium Imaging Survey (MIS) mode, which 
reaches a depth of 22.6 and 22.7 AB in FUV and NUV re- 
spectively, GALEX detects more than 90 percent of bright 
(r < 16.8), nearby (0 < z < 0.11) early-type galaxies ob- 
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served by the SDSS (see Figure 1 in Kaviraj et al. 2006). 
Given the large-scale nature of the GALEX dataset and the 
excellent overlap with existing optical surveys, it is there- 
fore possible to incorporate UV photometry into age-dating 
analyses of stellar populations. 

In this study we combine GALEX broad-band UV pho- 
tometry in two passbands (NUV: ~ 2310 angstroms and 
FUV: ~ 1530 angstroms) with optical photometry in the 
standard U, B, V, R, I filters to simultaneously estimate the 
ages and metallicities of GCs in M31. The GCs are drawn 
from the catalogue of Rey et al. (2005) who have combined 
GALEX imaging of th e M31 system with the optical cata- 
logue of lBarmbv et all J2000). The use of GCs is motivated 
by the fact that their star formation histories (SFHs) can 
be adequately parametrised by simple stellar populations 
(SSPs) of a given age and metallicity alone. They are there- 
fore ideal test objects to demonstrate the usefulness of UV 
photometry in constraining the age-metallicity parameter 
space. Although we derive both age and metallicity, our em- 
phasis is on the quality of the age estimation, as this is of 
primary importance in the determination of formation histo- 
ries. Note that, while the stellar models used here (described 
below) assume solar abundance ratios, GC age estimates 
derived using optical colours tend to be lower by ~ 8 per- 
cent if alpha-enh anced ([a/Fe] = 0.3) ratios are considered 
l|Kim et al.ll2002h . 

In the following sections, we compare the quality of the 
constraints on age and metallicity with and without UV 
photometry, investigate how the constraints change as vari- 
ous filters are excluded from the parameter estimation and 
compare our purely photometric analysis to previous stud- 
ies that have used spectroscopic indicators to estimate SSP- 
equivalent ages and metallicities for a variety of systems. 

2 PARAMETER ESTIMATION 

We begin by outlining the method used for parameter esti- 
mation. For a vector X denoting parameters in the model 
and a vector D denoting the measured observables (data) , 



prob(X|D) oc prob(_D|X) x prob(X), 



(1) 



where prob(X|D) is the probability of the model given 
the data (which is the object of interest), prob(D|X) is the 
probability of the data given the model and prob(X) is the 
prior probability distribution of the model parameters. As- 
suming a uniform prior in our model parameters and gaus- 
sian errors gives 

prob(_D|X) oc exp(-x 2 / 2 ), ( 2 ) 

where exp(— x 2 /2) is the likelihood function and \ 2 is 
defined in the standard way. prob(X|D) is a joint PDF, 
dependent on all the model parameters. To isolate the effect 
of a single parameter X\ in, for example, a two parameter 
model (prob(X|D) = prob(Ai , X2 \D)) we can integrate out 
the effect of X2 to obtain the marginalised PDF for Xi : 



prob(Ai|L>) = / prob{X 1 ,X 2 \D)dX 2 . (3) 
Jo 

In our analysis, the parameters X in the model are 



the age (t) and metallicity (Z) of the SSP in ques- 
tion, and the observables are colours constructed from 
FUV, NUV, U,B,V,R,I photometry. We choose a uniform 
prior in t, in the range to 14 Gyrs. We set the maxi- 
mum allowed age to 14 Gyrs using the WMAP estimate 
for the age of the universe (13.7 Gyrs, Spergel et al. 2003) 
- we note that structure formation probably commences 
~ 0.2-0.5 Gyrs after the creation of the Universe. For the 
metallicity Z, we use a uniform prior in the range -2.3 to 
0.36 dex, which are essentially the metallicity limits of the 
stellar models. For each GC analysed, we marginalise the 
age and metallicity parameters, take the peak (xo) of the 
marginal PDF P(x) as the best estimate and define one- 
sigma limits x+ and x_ as P(x)dx = 0.34 P(x)dx 
and f*° P(x)dx = 0.34 f*'^ P(x)dx respectively. 

The uncertainty in the parameter is then (x+ — X-). 
Note that we do not consider correlated errors in the two 
quantities. 



3 AGE AND METALLICITY CONSTRAINTS 

Young stellar populations (< 1 Gyr old) emit a substan- 
tial proportion of their flux in the ultraviolet, making UV 
photometry a good indicator of their presence. However, 
UV flux can also be generated by old, evolved stages of 
stellar evolution. Due to their longevity, hot, core Helium- 
burning horizontal-branch (HB) stars and their progeny are 
thought to be the dominant UV sources in old stellar popu- 
lations, such as GCs and elliptical galaxies. While metallicity 
is the dominant parameter controlling the HB morphology 
of Galactic halo-like stellar populations, the second parame- 
ter governing the pr operties of the HB has been the subject 
of mu ch debate (e.g. ISaraiedini et al.| [l997; Bcll azzini et al.l 
2001), concensus favouring age differences between clusters 
as t he major second parameter affecting HB morphology 
(e.g.|l]3|l994|). 

The stellar models used in this study employ detailed 
prescriptions for the systematic variation of HB morphology 
with metallicity and age, and are well-calibrated to Galactic 
GC populations and the integrated UV-op tical photome- 
try of local elliptical galaxies (see |yJ 120031 . and references 
therein). Main-sequence to Red Giant Br anch isochrones 
employed in these models are described in lYi et alj (1200 if) 
and the HB stellar tracks are des cribed in Yi et alj (1 19971 ). 
The flux library employed is from iLeieune et aL l | 19981 ). 

In Figure [TJ we illustrate the onset of UV flux in the 
models with the appearance of horizontal branch stars, as 
the stellar population evolves through intermediate ages of 

4 and 7 Gyrs to an age comparable to stars in the Galac- 
tic halo (~ 13 Gyrs). The age estimation procedure essen- 
tially fits the observed shape of the t/V-optical continuum 
to the models. In this context we note that, although the 
shape of the UV spectrum from a young and old population 
can appear similar, the slope of their optical continua are 
very different, which allows them to be differentiated easily. 
Thus, with access to a reasonably large set of UV and op- 
tical filters, age determinations of stellar populations such 
as GCs can be performed with reasonably good accuracy, as 
we demonstrate below. 

The purpose of this Letter is not to perform a rigor- 
ous study of the entire M31 GC system but to demonstrate 
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Figure 1. UV flux from old populations due to the onset of the 
horizontal branch population (top to bottom). Intermediate age 
populations (~ 4-7 Gyrs old) output negligible amounts of UV 
flux. The contributions from the MS-RGB and HB phases are 
shown in red and blue respectively. In all panels [Fe/H] = -1.7 
dex. The HB develops at higher (lower) ages for [Fe/H] greater 
(lower) than -1.7 dex. For example, for [Fe/H]~-1.2 and [Fc/H]~- 
2.1 dex, the HB develops at 12 and 8 Gyrs respectively. Note that 
this implies that metal-rich GCs will have weaker age constraints 
due to the later rise of the HB. 



the usefulness of UV photometry in reducing the AMD and 
present a method that exploits this property. We therefore 
extract parts of Rey et al. (2005) 's catalogue that best sat- 
isfy this aim - GCs with good detections in GALEX, which 
have the full set of optical U, B , V, R, I photometry and lie in 
regions of relatively low Galactic reddening (Eb-v < 0.2). 
We find 42 objects matching this description in the cata- 
logue. 

We begin by investigating the age and metallicity con- 
straints achieved using seven broad-band colours: FUV — V , 
NUV-V, FUV-NUV, U-B, B-V, V-R, V-I (Figure 
[2]) and compare them to the scenario where the UV colours 
are absent (Figure[3)l . The large AMD due to using only opti- 
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Figure 2. Age and metallicity constraints from UV and opti- 
cal photometry for two GCs in our sample. LEFT COLUMN 
Best-fit position and one-sigma contour. MIDDLE COLUMN 
Marginalised metallicity PDF (normalised). RIGHT COLUMN 
Marginalised age PDF (normalised). Top row shows a young clus- 
ter (148-200) and bottom row shows an old cluster (166-000). The 
gree n region shows a spect ros copic metallicity mea surement taken 
from lBarmbv et al.l j200Ch or lPerrett et al.1 hoO$\ . 
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Figure 3. Same as Figure [2] but using optical photometry only. 



cal passbands is clearly visible in Figure [3] The effect on the 
age estimation is severe - since the parameter space inside 
the contour is indistinguishable in terms of the goodness of 
fit (the best-fit value has very limited significance), the age 
of either GC can be anywhere between 2 and 13 Gyrs. The 
age constraint is far superior when UV photometry is added. 
The young cluster (148-200) now has a well-constrained age 
of 0.5-2 Gyrs and the old cluster has an age of 10-13 Gyrs. 
The photometry suggests that both clusters are metal-poor 
([Fe/H] < —0.5), in agreement with the spectroscopically 
determined metallicity (green region). 

In Figure [4] we summarise the age and metallicity con- 
straints from the full set of UV and optical photometry for 
the 42 GCs in our sample. The age constraints are robust, 
with the full extent of the age uncertainties lying between 
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Figure 4. Age and metallicity constraints from the full set of UV 
and optical photometry from our 42 GCs. 
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tometry may not always be available. For example, 
quality photometry is far more difficult to achieve in the U 
band than in longer optical wavelengths. It is therefore in- 
structive to compare the quality of the age and metallicity 
constraints when a subset of these passbands is excluded 
from the parameter estimation. We quantify this compari- 
son in Table 1. Since using only optical filters always gives 
the largest uncertainties, we compute the ratio of the uncer- 
tainty from using a particular set of filters to the uncertainty 
when only optical filters are used. We show the average val- 
ues of these ratios, for the age constraints (ft - column 9) and 
for the metallicity constraints (fz - column 10). A smaller 
value thus indicates a better constraint and smaller uncer- 
tainty in the estimation of age or metallicity. The filter sets 
are listed in order of the quality of the age constraint i.e. the 
first row indicates the filter set with the best age constraint. 
The constraints are virtually identical as long as either the 
R or I filters are present - we have therefore put them in 
the same column in Table 1. 

Looking first at constraints on the age estimates, we find 
that, for the full set of FUV, NUV, U, B, V, R, I filters (row 
1), the size of the marginalised age errors are ~ 10 percent of 
the optical-only case (f t = 0.12) i.e. the age estimates are ~ 
90 percent better constrained than those from using optical 
filters alone (row 7). Removing the longer wavelength R,I 
bands (row 2) or the shorter wavelength U band (row 3) 
slightly weakens the age constraint compared to the full set 
of filters. If only one UV filter is available, it is preferable to 
combine the FUV passband with optical photometry (row 
5) than the NUV passband (row 6). 

Similarly, the comparison of metallicity estimates (val- 
ues of fz in Table 1) indicates that adding UV photome- 
try to the optical U, B, V, R, I filters (or a subset of them) 
improves the metallicity estimates by approximately ~ 60 
percent compared to the optical-only case, fz refers to the 
underlying met allicity Z, conve rted from the [Fe/H] value 
using Table 2 of lYi et all I^OPl! ). 



Figure 5. TOP PANEL: Comparison of the best-fit ages of 
GCs to their (NUV - V)o and (FUV - V)o colours. BOTTOM 
PANEL: Comparison of the best-fit ages of GCs to their (U — B) 
and (B — V) colours. 



0.1 and 3 Gyrs. The metallicity constraint is not as robust 
as the age uncertainties, which is expected since the UV is 
a better indicator of age than metallicity. We find that, al- 
though the photometric metallicity constraint overlaps with 
its spectroscopic equivalent in all cases, the spectroscopic 
constraint is invariably more robust. 

In Figure [5] we compare the best-fit ages of each GC 
to their photometric colours. Contrary to the optical sce- 
nario, we can use the best-fit ages in this case because 
the age uncertainties are small. We find that GCs with 
well-constrained young ages tend to be bluer in both the 
(NUV — V)o and (FUV — V)o colours. In contrast, the com- 
parison to optical colours is more ambiguous - while there 
is an indication that the young GCs are bluer in (B — V)o, 
there is no clear separation between the young and old pop- 
ulations, as a result of the AMD. 

In practice, a full set of FUV, NUV, U, B, V, R, I pho- 



4 COMPARISON WITH SPECTROSCOPIC 
METHODS 

We now compare the quality of our photometric age 
and metallicity constraints with three studies which de- 
rive these pa r amete rs using spectroscopic line indices. The 
iTrager et all i|2000T ) study uses the H , Mg b, Fe5270 and 
Fe5335 indices to obtain tightly-constrained SSP ages for a 
sample of local elliptical galaxies, with age uncertainties of 
~ 0.5 — 1.5 Gyrs for young objects and ~ 2 — 3 Gyrs for old 
and inter mediate age objects ( see th eir Table 6A) . In a sim- 
ilar work, iProctor fc Sanson] l|2002h combine line strength 
measurements of 24 age and metallicity sensitive indices for 
the central regions of 32 galaxies in both the Virgo cluster 
and the field. Their age uncertainties are in the range ~ 1 — 3 
Gyrs (see their Tab l e 10 or Figure 12). In a more recent 
work, iBeaslev et ail (I2004T) measure 24 indices in the Lick 
system l|Worthev fc Ottavianil [l997l ; ITrager et alj|l998h for 
a sample of 30 globular clusters in M31. The uncertainties 
in their age-estimates are ~ 0.2-1.25 Gyrs (see their Figure 
12). Looking at Figure [3] we see that the uncertainties in 
our photometric age-determinations are ~ 0.1 — 2 Gyrs for 
the young clusters and ~ 0.5 — 3 Gyrs for the old clusters. 
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Table 1. Comparison of the quality of age and metallicity constraints for different filter sets, ft (fz) is the ratio of the age (metallicity) 
uncertainty using the filters in a given row to the age uncertainty using optical filters only (row 7). The filter sets are listed in order of 
the quality of the age constraint i.e. the first row indicates the filter set with the best age constraint. 
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We therefore find that our photometric analysis produces 
comparable age constraints to the spectroscopic studies. 

However, our metallicity constraints are arguably worse 
than those in the spectroscopic studies mentioned above. 
While we do not resolve the metallicity to better than 
±0.3 to 0.5 dex for any cluster, spectroscopic constraints 
on metallicity are generally tighter, of the order of ~ 0.1 
dex. This is expected, not only since the UV passbands are 
better indicators of age than of metallicity but because the 
combination of multiple metallicity-sensitive indices in these 
studies results in better metallicity constraints than can be 
achieved with photometry alone. 



5 CONCLUSIONS 

In this study we have demonstrated that the combina- 
tion of FUV (~ 1530 angstroms) and NUV (~ 2310 
angstroms) photometry with optical observations in the 
standard U, B, V, R, I filters can efficiently break the AMD 
caused by using optical filters alone. We have estimated ages, 
metallicities and their associated uncertainties for a sample 
of 42 GCs in M31 for which we have the full set of UV 
and optical photometry. We have found that using the full 
set of FUV, NUV, U, B, V, R, I filters produces age estimates 
that are ~ 90 percent more constrained and metallicity es- 
timates that are ~ 60 percent more constrained than those 
produced by using optical filters alone. The quality of the 
age constraints is comparable to those achieved using a large 
number of spectroscopic line strength indices. We find that if 
only one UV passband is available, then combining the FUV 
rather than the NUV with optical photometry is preferable 
in terms of constraining ages. Finally, while there is overlap 
between our photometrically determined metallicity ranges 
and the corresponding constraints from spectroscopic stud- 
ies for GCs in our sample, the photometric metallicity con- 
straints are generally weaker than their spectroscopic coun- 
terparts. 

We note that, since the UV is produced by either young 
(< 2 Gyr old) or evolved (> 10 Gyr old) stars, it is insen- 
sitive to intermediate-age populations, leading to a lack of 
intermediate-age GCs in Figure 4. To conclude, we suggest 
that broad-band UV photometry from the GALEX mission 
(which, at completion, will survey the entire local universe) 
could be used as an effective age indicator to break the 
AMD, and as a useful proxy for spectroscopic methods in 
constraining the age-metallicity parameter space. 
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